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Electroweak penguins and SUSY K^-K^ mixing with Neuberger quarks*^ 
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We present results for A7=3/2 and A5'=2 matrix elements relevant for CP violation in K mr decays and 
for the Ks—Kl mass difference in the standard model and beyond. They were obtained with Neuberger fermions 
on quenched gauge configurations generated with the Wilson plaquette action at /3=6.0 on an 18^ x 64 lattice. 



1. Introduction 

CP violation in K ^ irn decays and the Ks- 
Kl mass difference provide important constraints 
on the standard model (SM) and extensions, once 
the relevant hadronic matrix elements are accu- 
rately computed. In the following we present re- 
sults for the chirally leading A/=3/2 contribu- 
tions to direct CP violation, and for the AS'=2 
matrix elements which contribute to K^-K^ mix- 
ing in SM extensions. These are obtained using 
the same Neuberger propagators as for the spec- 
troscopy work presented in [1], where details of 
the simulation are given. We work with degen- 
erate u, d and s quarks of bare masses aruq — 
0.03, 0.04, 0.06, 0.08, 0.10. The chiral symmetry 
of Neuberger fermions guarantees that the mixing 
on the lattice is identical to that in the continuum 
and that our results are fully 0(a)-improved. 

2. Electroweak penguins 

At leading chiral order, the A/=3/2 contribu- 
tion to direct CP violation in if — > tttt is deter- 
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mined by the matrix elements ((7r7r)/=2|Q7,8|-f^''): 
of the electroweak penguin operators Qy.s- Soft 
pion theorems can be used to reduce one of the 
pions, so that ^ 



(1) 



{{nn)i=2\Q7,8\K°) (X ±{n+\Q'/^\K+) 

in the chiral limit, where is the chiral limit 
value of the pion decay constant {Ft^ — 92 MeV) 
and where Q^^g are the A/=3/2 components of 
Q7,8, given by 



.3/2 
!7,8 



- [{sd)v-Aiuu)v+A + {su)v-Aiud)v+A 



-{sd)v-A{dd)v+A] 



(2) 



with Q^^^ color diagonal and Q'^'^ color mixed. 
Thus, we calculate {7r^\Qj^^\K^) on the lattice 
using Neuberger fermions at finite quark mass 
and extrapolate the results to the chiral limit. 

Because we work in the isospin limit, eye 
contractions cancel and power-divergent mixing 
with lower dimensional operators is absent. This 
greatly simplifies the calculation. We begin by 
constructing the following ratios of correlation 
functions (a = 1): 



3/2 



7,8 



7,8 



(3) 



^The proportionality constant is convention dependent as 
is the overall phase of {-iT^\Q^^a\K^). 
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Figure 1. Dependence on initial time of the 
correlation function ratio Bpfr ^^r a variety of 
sources and sinks at fixed final time tf = 16 and 



for aniq = 0.04. 



with N7=3 and Ns=l. B^^^ measures deviation 
from the VSA of {n+\Q^^(^\K+) in the chiral limit. 
In Eq. (3), J''"^ = q'Tq with q = {I - aD/2p)q. 
Similarly, Q^^g are given by Eq. (2), with quark 
fields q replaced by q and with the appropriate 
rewriting in terms of Euclidean Dirac matrices. 

To quantify possible unwanted contributions 
from finite- volume zero modes, we vary the source 
and sink by letting Tij be one of 75, 7075, (1±75) 
or 7o(l ± 75). As shown in Fig. 1 for our next to 
lightest quark mass, amq = 0.04, which corre- 
sponds roughly to nis/i, we see no dependence 
of 'Bp^p. on source and sink in the fit region. 
Sp^p. displays very similar behavior. We inter- 
pret this as meaning that zero-mode effects are 
negligible around the kaon mass. The fit region, 
43 < < 51 and 13 < < 21, is chosen such 
that the ratios of Eq. (3) are asymptotic and free 
from wrap-around contributions. Sj^'^^^ yield the 

results for B^^^'''"'^'^ with the smallest error bars 
and we take these as the starting point for subse- 
quent analysis. 

Before proceeding with the chiral extrapola- 
tion of the matrix elements, we have chosen to 
renormalize them. We do so non-perturbatively 
in the RI/MOM scheme following [2]. Q7 and 
Qs mix under renormalization. We fix gluon 
configurations to Landau gauge and numerically 
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Figure 2. p^-dependence oi TZf^ {0 , p'^ ) . 

compute the relevant amputated forward quark q 
Green functions with legs of momentum p ~ \f]?". 
Aqj g{'mq,p'^). Then we determine the renormal- 
ization constants by requiring that the renormal- 
ized vertex functions have their tree-level values. 
Thus, we define the ratio: 



Tr {AvifUq, p^)Vvy 



(4) 



Tr |AQ^(m,,p^)PQj 

where i,j € {7,8} and the Vo are normalized 
projectors onto the spin-color structure of tree- 
level O = Q7,s,V- We extrapolate these ratios 
to niq = and fit the results to the OPE form 
(including discretization error terms) [3] : 



A 



ik 



kj 



+B,,{apf 



(5) 



where U^{p'^) describes the running of the renor- 
malization constants in the RI/MOM scheme, 
Zf^{p'^) = Z^^'Ufj^{p^), implemented at 2-loops 
[4,5] with Nf = Q and from [6]. For B- 
parameters, of course, the VSA's must be appro- 
priately renormalized. 

In Fig. 2 we show the p^-dependence of the ra- 
tios 7egi(0,p2) and the fits, for 1.5 GeV < p < 
3.2 GeV, to the OPE expressions of Eq. (5) to the 
orders displayed in the caption. 

We now turn to the chiral extrapolation. Here, 
complications arise due to the fact that we are 
working in the quenched approximation. Indeed, 
the VSA's of {'K+\Qy'^\K+) are ill-defined in the 
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gure 3. {Tr+\Q-'/g\K+)/F^ in the NDR scheme 
2 GeV are chiraUy extrapolated. 



{K^\0^\K°) with 

Oi = [sd]v-A[sd]v-A 
O2.3 = [s(i]s-p[sd]s-p 
04,5 = [sc?]s-p[srf]s+p 

where: 



(unmix, mix) (9) 
(unmix, mix) 



• {K"\Oi\K°) is the SM contribution, best 
given in terms of 



B 



K 



3 (ifoiOili^") 



16 



chiral limit due to the presence of quenched chi- 
ral logs [7]. We could consider {■K+\Qy^l\K+) / F'^ , 
with F the decay constant of our degenerate 
and 7r+. Instead, we choose to extrapolate 
{■K^\Qy^\K^) / F'^ , which have similar mass de- 
pendences in the quenched and unquenched the- 
ories. Indeed, at 1-loop in xPT, 



K+) 



for Nf^3 [8,9] andNf^Q [7,9]. This similarity in 
mass dependence may be an indication that some 
quenching errors cancel in these ratios. These ra- 
tios have the added advantage that they are free 
of chiral logarithms at 1-loop, and should hence 
have smoother chiral behaviors. We show the cor- 
responding chiral extrapolations in Fig. 3, where 
we have added to Eq. (6) terms proportional to 
(M/47rf-^)^ to account for higher orders. We ob- 
tain, in the NDR scheme at 2 GeV 

^ = 1.9±0.3±??GeV^ 



lim 



lim 



r+\QT\K^ 
F2 



(7) 



= 12. ± 2. ± ?? GeV^ , (8) 



where ?? stands for systematic errors which have 
yet to be determined. We postpone to a later pub- 
lication the comparison of our results with non- 
lattice [10-13] and quenched Wilson [14], domain- 
wall [15,16] and overlap [17] lattice results. 

3. AS=2 transitions beyond the SM 

In extensions of the SM, analysis of K'^- 
K'^ mixing generically requires knowledge of 



. {K^\05a\K^) 
rriu = nrid 



2y^{-n+\Qy^l\K+) for m. 



• {K^\Oi\K^) is 0{p^) in the chiral expansion 
for i ^1. 

The last point indicates that non-SM matrix ele- 
ments are expected to be larger than {K^\Oi \K^), 
which is 0{p^). To make this explicit, we define 
the ratios {i=2, • • • , 5) 



r F'^ 1 




{K°\0^\K°)' 




L ^"kJ 


expt 


.F2 {K0\Oi\K0)_ 


lat 



,(10) 



where BSM stands for "beyond the SM" and 
where M and F are the mass and "decay con- 
stant" of the lattice K'^. These ratios have a 
number of advantages: 



Hat 



Hat 



is dimensionless 

is finite in the chiral limit 



• Rf^^{M'^) measures directly the ratio of 
BSM to SM contributions. 

We obtain the required matrix elements from 
ratios of 3-point to two 2-point functions, as we 
did for (TT+IQy^g Pseudoscalar sources and 

sinks were used, except for {K°\Oi\K^) where 
left-handed current sources and sinks were taken 
to completely eliminate zero-mode contamina- 
tion [3]. We also perform an RI/MOM non- 
perturbative renormalization as above: Oi renor- 
malizes multiplicatively; 02,3 mix; the O4 5 mix- 
ing matrix was already shown. 

In Fig. 4, we present new results for 
B^^^{2GeV) as a function of pseudoscalar mass 
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Figure 4. Bf^^{2GeV) vs at /3 
18^ X 64 and 16^ x 32 [3] lattices. 
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Figure 5. Rf^'^{A'P) in the RI/MOM scheme at 
2GeV vs JVP. The vertical line marks M'^=Mj.. 



squared M^, together with our published results 
obtained at the same lattice spacing but on a 
smaller 16^x32 lattice [3] . The results are entirely 
compatible, indicating that finite-volume elfects 
are small around the physical kaon mass. We find 
Bf^^{2GeV) = 0.60(5)(1), which is to be com- 
pared with the value of 0.63(6)(1) that we found 
on the smaller lattice. Our results are compatible 
with the quenched benchmark result of JLQCD 
[18], as well as with the overlap results of [17]. 

In Fig. 5 we plot results for Rf^^iM"^) in the 
RI/MOM scheme at 2 GeV as a function of M^. 
We find that the BSM matrix elements are en- 
hanced at M—Mk by factors ranging from ap- 



proximately 5 to 20 compared to the chirally 
suppressed SM one. This enhancement is sig- 
nificantly larger than the one observed in the 
quenched, Wilson fermion calculation of [14]. 

4. Outlook 

We are investigating systematic errors on the 
results presented above. In particular, we are per- 
forming the same calculation at /3 = 5.85 on a 
14'^ X 48 lattice to quantify discretization errors. 
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